Draft version March 2, 2013 

Preprint typeset using I^TeX style emulateapj v. 5/2/11 



CORONAL MASS EJECTION MASS, ENERGY, AND FORCE ESTIMATES USING STEREO 
Eoin P. Carley 1 , R. T. James McAteer 2 , and Peter T. Gallagher 1 

Draft version March 2, 2013 

ABSTRACT 

Understanding coronal mass ejection (CME) energetics and dynamics has been a long-standing 
problem, and although previous observational estimates have been made, such studies have been 
hindered by large uncertainties in CME mass. Here, the two vantage points of the Solar Terrestrial 
Relations Observatory (STEREO) COR1 and COR2 coronagraphs were used to accurately estimate 
the mass of the 2008 December 12 CME. Acceleration estimates derived from the position of the CME 
front in 3-D were combined with the mass estimates to calculate the magnitude of the kinetic energy 
and driving force at different stages of the CME evolution. The CME asymptotically approaches a 
mass of 3.4 ± 1.0 x 10 15 g beyond ~10 R@. The kinetic energy shows an initial rise towards 6.3 ± 3.7 x 
10 29 erg at ~3i?©, beyond which it rises steadily to 4.2 ± 2.5 x 10 30 erg at ~18i?©. The dynamics 
are described by an early phase of strong acceleration, dominated by a force of peak magnitude of 
3.4 ± 2.2 x 10 14 N at -3 i? , after which a force of 3.8 ± 5.4 x 10 13 N takes affect between -7-18 R Q . 
These results are consistent with magnetic (Lorentz) forces acting at heliocentric distances of <7i?©, 
while solar wind drag forces dominate at larger distances (>7 R®). 
Subject headings: Sun: corona - Sun: coronal mass ejections (CMEs) 



1. INTRODUCTION 

Despite many years of study, the origin of the forces 
that drive coronal mass ejections (CMEs) in the solar 
corona and interplanetary space are not well understood. 
From an observational viewpoint a complete understand- 
ing of CME kinematics, dynamics and forces requires 
not only a study of CME speed, acceleration and expan- 
sion but also an accurate knowledge of CME mass. The 
measurements of CME mass combined with acceleration 
measurements can be used to quantify the magnitude of 
the force that drives a CME. Knowledge of this force 
magnitude can lead to an identification of the possible 
origin of the CME driver. 

There are numerous theoretical models that attempt 
to explain the triggering of CME eruption and its con- 
sequent propagation. Each describe the destabilization 
and propagation of a complex magnetic structure, such 
as a flux rope , via mechanisms that include the catas- 



trophe model (Forbes fc Isenberg 1991 [Forbes fc Priest 
tiochos et al. 1999; Lynch et al. 2008) , or a toroid al in- 



1995 ; [Lin fc F orbes 2000), magneticTJfeakout model (|An- 



stability model ([Chen 1996; K liein^Toro^[2006| ) . The 
loss of equilibrium induced by such mechanisms results in 
CME propagation into interplanetary space. The predic- 
tions of these models have been investigated in observa- 
tional studies whereby the CME kinematics are used to 
constrain what forces might be at play and hence which 
model best describes CME propagation. Such studies 
show that early phase propagation can be reasonably de- 
scribed by the existing models (or a combinatio n of them) 
involving some form of magnetic CME driver flManoha- 



ran fc Kundu| [20031 |Chen et al.|[2006| |Schrijver et"aT 



2008; Lin et al.|2010 ), and that aerodynamic drag of the 



solar wind may have a significant role at later stages of 
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CME propagation ([Howard et al."||2007| |Maloney fc Gal- 



lagher 



2010 Byrne~t al.||2010|). Comparisons between 



modeling and observational estimates of the forces that 
drive CMEs requires an accurate determination of CME 
kinematics properties as well as CME mass. 

To date, the most prevalent method of determining 
CME mass has been through the use of white light 
coronagraph imagers, such as the Large Angle Spec tro- 
scopic Coronagraph (LASCO; |Brueckner et al.||l995[) on 
board the Solar and Heliospheric Observatory [SO HO] 
Domingo et al.|1995[ ) and the twin Sun Earth Connection 



Coronal and Heliospheric Investigation (SEC C HI) COR1 
and COR2 coronagraphs ( jHoward et al.|[2008| ) on board 
the Solar Terrestr ial Relations Observatory [STEREO] 
Kaiser et al.||2008 ). The white-light emission imaged by 



such coronagraphs occurs via Thoms on scattering of pho- 
tospheric light by coronal elec trons ( |Minnaert| [l930 van 
de Hulst 1950 Billings| p"966), the so called K-corona. 
From classical Thomson scattering theory, the intensity 
of the light detected by an observer depends on the parti- 
cle density of the scattering plasma. Hence, any density 
enhancement, such as a CME, over the background coro- 
nal density appears as enhanced emission in white light. 
The enhanced emission allows for a calculation of the 
total electron content and hence mass. 

Some of the first measurement of CME mass using scat- 
tering theory were c arried out by Munro et al. ( 1979 ) and 



Poland et al. ( 1981 ) using space-based white light corona- 
graphs on board Sky Lab and U.S. military satellite PI 8-1. 
Both the early studies and later statistical investigations 
determined that the majo rity of CMEs have masses in 



the range of 10 13 -10 16 g, dVourlidas et al. 2002 2010) 



However, due to only a single viewpoint of observation, 
the longitudinal angle at which the CME propagates out- 
wards was largely unknown in these studies and it is gen- 
erally assumed that the CME propagates perpendicular 
to the observers line-of-sight (LOS). There is also the 
added assumption that all CME mass lies in the two- 



2 



dimensional plane-of-sky (POS). Such assumptions can 
l ead to a mass u nderestimation of up to 50% or more 
(Vourlidas et al. 2000). More recent studies have em- 
ployed the two viewpoint capabilities of the STEREO 
mission to determine t he mass of numerous CMEs w ith 
much less uncertainty ( |Colaninno fc Vourlid as 2009). 

In this paper, we analyze mass development of the 
2008 December 12 CME using the STEREO COR1 and 
COR2 coronagraphs.We use a well constrained angle of 
propagation to determine the mass and position of the 
CME. Combining the mass measurements with values 
for CME velocity and acceleration, the kinetic energy 
and the magnitude of the force influencing propagation 
is determined for each point in time. Section 2 describes 
the observations of the event from first appearance of the 
front in COR1 A and B to the time when the front exits 
the COR2 A and B fields of view. Section 3 describes the 
methods by which the mass, energy, and force are calcu- 
lated with a priori knowledge of the propagation angle. 
Section 4 includes the results and Section 5 discusses the 
possible forces attributable to the observed accelerations 
and whether they are magnetic or aerodynamic in origin. 
This is followed by conclusions in Section 6. 

2. OBSERVATIONS 

The COR1 images used in this analysis span from 
2008 December 12 04:05 UT to 15:45 UT, with a ca- 
dence of 10 minutes. The three polarization states of 
COR1 were combined to make total brightness images 
in units of mean solar brightness (MSB). Base difference 
images were produced using the 04:05 UT image (in both 
CORlAandB) as a background to be subtracted from 
all subsequent images. A sample of such images for both 
COR1 A and B can be found in Figure [TJ The COR2 
images analyzed range from from 07:22 UT to 17:52 UT, 
with a cadence of 30 minutes. As with the COR1 images, 
total brightness images were created for COR2, and a set 
of base difference images were then produced using the 
07:22 UT image as a suitable background. A selected set 
of images from COR2 can be found in Figure [2] 

At 04:35 UT the leading edge of a CME appeared in 
COR1 A and B coronagraphs at a height of ~1.4i2 , off 
the east and west limb respectively. In COR1 B the CME 
first appears as a set of rising loop-like structures followed 
by a prominence, part of which appears to fall back to 
the surface at 08:00 UT while the remainder was ejected 
and follows the rising loop-like structures which eventu- 
ally become the CME front. The rising prominence was 
not apparent at any stage of the propagation in COR1 A 
and the advancing front remains the only distinguishable 
facet of the CME from this line-of-sight (LOS). 

A noteworthy caveat of using base difference imaging is 
the assumption that the background corona in the pre- 
event image has the same brightness in all subsequent 
images. This may not always be true and any excess 
brightness in the pre-event image will produce negative 
pixel values in the base difference. This is apparent in the 
COR1 images as the CME interacts with a streamer, dis- 
placing it as the leading CME front expands laterally as 
well as moves outward. The streamer is visible as a dark 
feature that grows with time at the southward flank of 
the CME in the COR1 B images, Figure [lj The black ar- 
eas are indicative of negative pixel values. The COR1 A 
images also suffer from negative pixels, especially at later 



times, see Figure [T] top row, 09:15 UT image. The front 
of the CME starts to exit both the A and B field of view 
at -08:35 UT. 

The CME first appears in the COR2 field of view at 
-07:52 UT with the CME apex at a height of ~3i? in 
both A and B images. In the B coronagraph, by 10:52 UT 
the three part structure of core, cavity, and bright front 
is clearly visible and the overall structure grows in size 
as the CME propagates to larger heights. The core be- 
comes more tenuous and the mass distribution becomes 
homogenous after 15:52 UT when the front starts to exit 
the field of view. The distinction between core and front 
is not as clear in COR2 A and the mass distribution ap- 
pears more homogenous throughout the propagation. As 
with the COR1 images, COR2A is also affected by ex- 
cess brightness in the pre-event image, as is apparent by 
a growing dark feature in its southern half. As the pre- 
event image for COR2 B is the cleanest of the pre-event 
images (it contains the least contamination by stream- 
ers), the COR2B data are considered the best candidate 
for accurate CME mass measurements. 



3. CME MASS MEASUREMENT METHODS 

The method by which mass measurements are derived 
from white light c oronagraph imag es is based on theory 
first developed by |Minnaert| (p. 930 ) in which the scatter- 
ing geometry of a single electron at a particular point 
in the solar atmosphere is considered. Fur ther develop- 
ment of the theory by van de Hulst| ( |1950| ) led to the 
derivation of what are now known as the van de Hulst 
coefficients. The coefficients treat each component of the 
incident electric field vector separately and take into ac- 



count the fin i te size of the solar disk (Minnaert 1930 



|Billings||1966| |Howard fc Tappin||2009| . An important 
tact arising from these expressions is the dependence of 
scattering intensity on the angle, \i between the radial 
vector from sun centre to the scattering electron and a 
position vector from observer towards the electron-the 
LOS, see Figure [3J Scattering efficiency is minimized 
when this angle is 90°. However, along the LOS such 
an angle occurs at the point of minimum distance from 
sun centre where the incident intensity (that the elec- 
tron receives) and electron density are maximized. This 
means scattered light in the corona is most intense along 
a plane perpendicular to the observer's LOS despite the 
efficien cy of scattering being mi nimized at such viewing 
angles ( jHoward fc Tappin||2009[ ) . This plane perpendic- 
ular to the LOS is known as the plane-of-sky (POS) 

Studies using single LOS coronagraph data are often 
hindered by the unknown CME propagation angle from 
the POS, e.g., unknown (or x) m Figure [3] This leads 
to the incorrect angle being used when inverting the van 
de Hulst coefficients to calculate the number of electrons 
contributing to the scattered light. Furthermore, because 
the 3-D extent of the CME is unknown it is also assumed 
that the CME is confined to the 2-D sky plan e, leading 
to a sign ificant CME mass underestimation (|Vourlidas 
etaT1[2000| . 
The CME 



IE of 2008 December 12 was Earth-directed 



(Byrne et al. 2010), making it roughly the same angu- 
lar distance from both the STEREO A andB spacecraft, 
then located ±45 degrees from Earth. This known angle 
of propagation was used to convert from pixel values of 




Fig. 1. — Selection of base difference images of the CME in COR1 A (top row) and COR1 B (bottom row), with pixel values of grams. 
The CME is quite faint in the A images and appears not to have as much structure as in B. There is a large contribution to mass from 
a near-saturated region to the upper flank of the CME in the B images. Such saturation in the mass images coincides spatially with the 
prominence in total brightness images. 



MSB to grams via the expression 



^pixel 



B P 



x 1.97 x 10" 24 g 



(1) 



where B b s is the observed MSB of the pixel, B e is the 
electron brightness calculated from the van de Hulst co- 
efficients, and 1.97xl0 _24 g is a factor that converts the 
number of electrons to mass, assuming a completely ion- 
ized corona with a composition of 90% hydrogen and 
10% helium. The known angle of propagation allowed 
the correct value of B e to be computed resulting in a sig- 
nificant reduction in the uncertainties associated with the 
propagation angle. The largest remaining uncertainty is 
the unknown angular width along the LOS. This uncer- 
tainty was quantified in a similar ap proach to the method 
outlined in |Vourlidas et ah] ( 2000 ) . This simulates the 
brightness of a CME with homogeneous density distribu- 
tion and finite angular width along the LOS-longitudinal 
angular width A9i ong , allowing calculation of a simulated 
observed mass. Comparing this to the actual mass al- 
lowed for an evaluation of CME mass underestimation 
for given values of A6i ong . Since the values for A0i O n a are 
unknown, the expression derived in |Byrne et al.| ( |2010[ ) 
for the latitudinal angular width of this CME as a func- 
tion of height, A0i at (r) = 25r 0,22 , was used to define an 
upper limit to A0i ong . It was assumed the CME longitu- 
dinal angular width is no more than twice the latitudi- 
nal angular width, or A0i ong ^2xA6i at . Such an upper 
limit is in agreement with simulations of flux-rope CMEs 



which give a typical aspec t ratio of broadside to a xial an- 
gular extents of 1.6-1.9 flKrall fc St. Cyr||2006| ). Hence 
the value for A0i ong at each height was used to obtain 
the simulated mass underestimation estimates described 
above. The heights and angular widths used in this study 
produced CME mass underestimation estimates of be- 
tween 5-10% for finite angular width uncertainty. An 
extra mass uncertainty of 6% was added to account for 
the assumption of coronal abundance of 90% hydrogen 
and 10% helium which can lead to slight er rors while 
converting from pixel values of MSB to grams ( |Vourlidas| 
et al.||2010|). 



To calculate the CME mass a user-selected area (the 
extent of the CME, for example) of the base difference 
image was chosen and the pixel values within this area 
were summed to obtain total mass. Figure [2] COR2 B 
images show an example of the sector over which pix- 
els were summed (the smaller sectors indicate a different 
summing region used at a later stage). The selected area 
was chosen for each image in the time sequence of CME 
propagation so as to determine the mass variation with 
height in COR1 and 2 using both A and B. The selection 
of an area by a point and click method is of course a sub- 
jective identification of the the extent of the CME, so it 
is susceptible to user- generated uncertainties. To quan- 
tify these uncertainties the mass was obtained for each 
coronagraph image in the time sequence (as described 
above) and the process was repeated five times in order 
to obtain the mean CME mass for each image and the 
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Fig. 2. — Selection of base difference images of the CME in COR2 A (top row) and COR2B (bottom row), with pixel values of grams. 
The CME is clearly distinguishable in both fields of view. Only the B field view shows clearly the three part structure of core, cavity and 
front. The COR2B images were used to measure core and front mass separately 

The deflection of a small streamer during CME prop- 
agation produces negative pixels in the base difference 
images. The effect is particularly apparent in the C0R1 
images, Figure [I] It is difficult to unambiguously distin- 
guish between streamer and CME, making it difficult to 
quantify the uncertainty introduced due to streamer in- 
teraction. To make an estimate of the streamer's effects, 
a calculation of its mass in the pre-event image was made. 
A number of different samples of the area of the streamer 
in the C0R1 B pre-event image that effects all subse- 
quent images produced a mass estimate of ~5xl0 14 g. 
This mass was used as a measure of the uncertainty in- 
troduced due to streamer interaction in the C0R1 B im- 
ages. A similar analysis of the C0R1 A pre-event images 
gave a streamer mass estimate of ~7xl0 14 g. C0R2 im- 
ages are relatively unaffected by significant changes in 
background coronal brightness and do not suffer from 
negative pixel values to as large an extent as C0R1. The 
pre-event image of C0R2 B is particularly clean and free 
of background streamers, hence C0R2 B images are con- 
sidered to provide most accurate CME mass estimation. 

Finally, in order to obtain a more complete and con- 
tinuous estimate of CME mass growth, the masses deter- 
mined from both C0R1 and C0R2 coronagraphs were 
summed in those cases where image times of the inner 
and outer coronagraphs overlappeqj The overlap in the 













c 


Sun \ 




POS 










. D 






[x 


LOS 



Fig. 3. — Schematic showing the relative orientation of the line- 
of-sight (LOS), and the plane-of-sky (POS). Electron position is 
at point P and C is Sun center. The vector CP may also repre- 
sent CME propagation direction. Scattering efficiency is heavily 
dependent on the angle (or x) an d is least efficient when = 0° 
(X = 90°). 

standard error on the mean. This standard error was de- 
fined as the uncertainty due to user bias in the point and 
click method of CME identification. The height at each 
measurement interval was taken to be the heliocentric 
distance of the CME apex in the image i.e., the apex of 
the front was chosen by simple point-and-click method. 
The uncertainty on the apex height was also found by 
the standard error on five runs. 



3 A difference in cadence of the inner and outer coronagraphs 
means that the images closest in time have a three minute off- 
set e.g., a COR1 image taken at 07:25 UT was considered to be 
coincident in time with the COR2 image at 07:22 UT 
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inner and outer corongraphs' fields of view was also taken 
into account in this summation. 

A concise measurement of the CME kinematics, such 
as velocity and acceleration , were taken from the results 
of the study of |Byrne et al.|pQ10[ ). Since these kinemat- 
ics take into account the true three dimensional surface 
of the front they provide reliable estimates of CME ve- 
locity and acceleration in 3-D space. These velocity and 
acceleration measurements were used in the calculation 
of kinetic energy and total force on the CME for each 
point in time. The CME mass used in all energy and 
force calculations was the asymptotic mass it approaches 
at later stages of its evolution beyond 10 Rq as observed 
from the STEREO B spacecraft i.e., 3.4 ± 1.0 x 10 15 g. As 
will be shown, there is good motivation for the use of con- 
stant mass in the magnitude of kinetic energy and force 
estimates. 

4. RESULTS 

4.1. CME Mass Estimates 

The results of the calculation for CME mass develop- 
ment with time and height for both STEREO A and 
B coronagraphs are shown in Figure [4j In panel (a), 
the height values are those taken from a point-and-click 
method of tracking the CME apex; these heights are cor- 
rected for CME propagation angle of ~45° . In both pan- 
els (a) and (b), the mass estimates of STEREO A and 
B follow a similar trend and have similar values at each 
stage in the propagation. Such good agreement between 
mass values is a good indicator that ~45° is the cor- 
rect angle of propagation from the sky-plane. A change 
in the cadence of mass measurements is noticeable at 
-08:00 UT (or >5i? ). This is due to the use of only 
COR1 images (with a cadence of 10 minutes) prior to this 
time, and the use of the COR2 plus COR1 images af- 
ter this time (the cadence of these measurements follows 
that of COR2- 30 minutes). Comparing A and B below 
4.5 i?0, mass values show a similar trend and increase at 
the same rate, but at approximately 3 Rq the mass mea- 
surements in COR1 B appear to increase to a much larger 
value then fall again. This effect is visible in the COR1 A 
measurements, albeit diminished. It is probably due to 
the presence of a prominence which contains a significant 
mass content and therefore contributes a large amount 
to total measured CME mass. Also, early on in its prop- 
agation, the prominence may still be emitting H-a line 
radiation (656.28 nm) due to the larger fraction of neutral 
hydrogen at its cooler temperatures. The COR1 imag- 
ing passband is centered on H-a so any emission in the 
prominence from neutral hydrogen could be contributing 
to light received by the COR1 coronagraphs, this is ap- 
parent from the saturation region in the COR1 B images 
in Figure [I] Since this is resonance line emission, and 
not Thomson scattered emission, it leads to an erroneous 
measurement in CME mass. Thus, it is assumed the 
larger rise and fall in CME mass is caused by the promi- 
nence entering and exiting the COR1 B field of view. The 
effect is diminished in COR1 A since the prominence does 
not enter the FOV to as large an extent as in COR1 B. 
The interpretation that the 'mass bump' is not actual 
mass growth (or loss) is supported by previous measure- 
ments where CME mass increase follows a trend with 
height described by M cme (h) = M (l - e~ h/ha ), where 
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Fig. 4. — CME mass development with height (a) and time (b), 
for the 2008 December 12 CME. After ~08:00 UT (>5i?©) the 
masses from the inner and outer coronagraphs are summed to 
show uninterrupted mass development from ~2-20 Rq over a pe- 
riod of 12 hours. The small bump in the CME mass at ~07:00UT 
(~4i?0) is probably due to an unknown amount of H-ck emission 
from the prominence. Mass of CME front and core are also shown, 
red 'X' and blue square, for COR2B, panel (b). After 14:52 UT 
they share approximately equal mass. The inset of (a) shows mass 
development with height for COR2B only; the red curve repre- 
sents a fit to the data whereby the mass asymptotically approaches 
3.4 ± 1.0 x 10 15 g. 

M a is the final mass the CME approaches asymptotically 
and h a is the height at which the CME reaches 0.63M a 



(Colaninno & Vourlidas 2009), with no 'bump' in mass 
earlier on. The decline in mass after the peak may be ex- 
plained by the ionization of neutral hydrogen such that 
H-a emission diminishes and simply becomes Thomson 
scattering of free electrons, as with the rest of the CME 
material. 

In order to produce a fit to the data, the COR2B 
mass results were chosen because its pre-event image 
was largely free of any bright streamers or other features 
which introduce unwanted effects in the production of 
base difference images, as described above. A fit with the 
above equation resulted in a final asymptotic CME mass 
of 3.4±1.0xl0 15 g, with a scale height of h a = 2.9 R Q . 
This fit is plotted along with the COR2 B data in the in- 
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set panel of Figure [4^ a). Note that the mass increase is 
due to material coming up from below the occulting disk, 
and not actual mass gain of the CME. The uncertainty 
on the above asymptotic mass value was taken to be 30%, 
from the largest uncertainty due to finite width, the con- 
version factor uncertainty as described above, the stan- 
dard error user-generated uncertainty, and uncertainty 
due to streamer interaction. 

In each image where the CME core and front are distin- 
guishable, their masses were measured separately. This 
was carried out by user selected regions demarcating the 
areas of core and front, see COR2B at 12:22 UT and 
14:52 UT in Figure [2] for an example of the separate core 
and front sectors over which pixel values were summed to 
obtain total mass. The uncertainties due to finite width 
of the observed object also apply to the core and front 
measurements, however, since the widths of these partic- 
ular areas of the CME are unknown we chose the max- 
imum uncertainty of 10% from the above analysis since 
neither core nor front can be any wider than the max- 
imum width assigned to this CME. The remaining un- 
certainties described above were also applied. The mass 
development of core and front with time is shown in Fig- 
ure gb). The two mass measurements are subject to an 
observational effect of apparent exponential mass growth, 
however by the time the CME is fully in the field of view 
at 14:45 UT the core and front share approximately equal 
mass. 

4.2. CME Forces and Energetics 

In the following calculations, all measurements of 
force and kinetic energy use the asymptotic mass of 
3.4 ± l.OxlO 15 g and not the instantaneous mass values 
calculated from each coronagraph image i.e., the CME is 
considered to begin its propagation with this mass and 
does not acquire any mass as it propagates. 

Estimates of the force and kinetic energy use the 3- 
D velocity and acc eleration measurements produced by 
Byrne et al. (2010). Their method firstly identifies the 
CME front in each coronagraph image using a multiscale 
edge detection filter. The front edges were then used to 
define a quadrilateral in space into which an ellipse is fit, 
this method is known as elliptical tie-pointing. This was 
done for multiple horizontal planes through the CME so 
that the fit ellipses outline a curved front in 3-D space. 
The speed and acceleration were then deduced from the 
change in position of the front, with time, through the 
STEREO COR1, COR2 and HI fields of view. Since 
mass measurements in this study use only the COR1 and 
COR2 coronagraphs, HI kinematics measurements have 
been excluded here. The CME front position uncertainty 
in STEREO A and B coronagraphs was determined from 
the filter width in the multiscale analysis. Velocity and 
acceleration uncertainties were then propagated from po- 
sition uncertainty. Figure [5|a) shows CME velocity as a 
function of heliocentric distance, along with acceleration 
in panel (b). 

The CME kinetic energy was calculated using E^in — 
l/2M cme ^ me , where M cme is the final asymptotic mass 
of 3.4 ± l.OxlO 15 g and v cme are the instantaneous veloc- 
ity measurements, results of this calculation are shown 
in Figure |5|c). The kinetic enerev shows an initial rise 
towards 6.3 ± 3.7x 10 29 ergs at ~3i?©, beyond which it 
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Fig. 5. — (a) CME velocity as a function of heliocentric distance, 
including a fit to the data produced u sing an aerodynamic drag 
model beyond ~7 R@ (| Byrne et al.|2010| >. (b) Acceleration of CME, 
including fit, derived from the velocity data and fit. Panel (c) and 
(d) show the kinetic energy and force, respectively, both calculated 
using constant CME mass of 3.4 ± 1.0 x 10 15 g and kinematics 
results from (a) and (b). Also shown are the fits to energy and 
force produced from fits to velocity and acceleration. 

rises steadily to 4.2 ± 2.5 x 10 30 ergs at ~18 these val- 



ues a re si milar to those reporte d in |Vourlidas et al. (2000 



2010]) and |Emslie et al . 
The total force on 



2004). 
CME 



total 



= M, 



cme^cme i 



the 
where 



was calculated using 
M cme is as above and 
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a cme is taken from the instantaneous acceleration val- 
ues. As shown in panel (d) of Figure [5| the force ini- 
tially grows significantly, reaching a maximum value of 
3.4±2.2xl0 14 N at ~3i? . The early rise and fall in 
acceleration (or force) is in agreement with a previous 
study of a CME observed to reach peak acceleration at 
~1.7i?0 after which it reaches a con stant velocity be- 
yond ^3.4 i? ( |Gallagher et al. 2003). Such results are 
also found in a statistical study which shows that the ma- 
jority of CMEs have peak acceleration in the low corona 
with a mean heigh t of maximum acceleration at 1.5 R® 



dBein et al.| 



Zhang et al. 



20111). Sim ilarly, observational studies by 
( |2001p and |Zhang et ah] fl2004D also show 
early phase peak acceleration between 2-5 Kq and forces 
on the order of 10 15 N and 10 12 N, depending on whether 
the CME shows large initial acceleration or a slow, more 
gradual acceleration. 

After this early peak, the force drops to an average 
value of 3.8±5.4xl0 13 N at distances between 7-18 Rq. 
It is apparent from Figure J5^a) that the velocity con- 
tinues to increase beyond 7 R@ > implying that a positive 
radial force must be present. To clarify this, a fit to the 
velocity data using a model for sol ar wind drag on th e 
CME beyond 7 R Q (as outlined in |Byrne et~aL| ( |2010[ )) 
is shown in Figure p|a). Although the data suggest a 
non-monotonic increase in velocity, the fit reveals that 
propagation is best described by a steadily increasing 
velocity between 7-18 Rq. The acceleration and kinetic 
energy curves derived from this velocity fit are shown in 
Figure |5|b) and (c) . In Figure |5|d) , the curve for the 
force derived from the velocity fit initially deviates from 
the data at ~7i?©, however beyond this distance there 
is good agreement with the data and the derived force is 
entirely positive. This suggests that the solar wind exerts 
a positive aerodynamic drag force on the CME, resulting 
in a velocity that approaches the asymptotic solar wind 
speed at large heliospheric distances. 

5. DISCUSSION 

It should be noted that Figure [4] shows an overall expo- 
nential increase in CME mass with height which could be 
interpreted as the CME rapidly gaining mass as it prop- 
agates. Care should be taken with this interpretation 
since this apparent exponential mass increase is almost 
certainly due to the CME moving into the field of view, 
therefore allowing us to measure more of its mass con- 
tent; such an inte rpretation is in agreem ent with similar 
assertions made in |Vourlidas et al. ( 2010 ) . It is difficult to 
distinguish between actual CME mass growth and an ap- 
parent growth due to more of the CME being observed. 
If the initial early rise in CME mass is assumed to be 
an observational artifact then we can interpret the CME 
mass to be in the range of (3-3.5) xlO 15 g for most of 
its early propagation i.e., the CME already has such a 
mass before launch and does not acquire more mass (via 
inflows or otherwise) during propagation. Such an inter- 
pretation is in agreement with CME mass measurements 
calculated from dimmings in STEREO Extreme Ultra- 
violet Image (EUVI) images, which show the mass cal- 
culated from EUV images to be approximately equal to 
CME mass in COR2 imag es, tueuvi /mcoR2 = 1.1 ±0.3 
flAschwanden et al.|[2009| Once the CME bubble is in 
the held oi view at ~1 Hq the mass in its entirety can 
be measured and the increase beyond this point, if any, 



is slow and steady, Figure [4j 

The early stages of CME propagation are dominated by 
a sharp rise to a peak force of 3.4 ± 2.2 x 10 14 N at ~3 R@ 
followe d by a sharp decline, Figure|5td) . The catastrophe 



orbes & Priest 



1995 



model dForbes fc Isenberg||1991 
Lin fc For bes 2000), magnetic 

chos et al. ||1999; Lynch et al. 2008ft, and toroidal insta- 



breakout model (fAntio 



bility model ( |Chen|19 96; Klien r^Torok|2006| ) employ a 
number of forces acting on the CME to produce an over 
all acceleration into interplanetary sp ace. For exa mple, 
the toroidal instability model used by ||Chen| p996) uses 
a Lorentz hoop force (or Lorentz self- force), solar wind 
drag, and gravity to provide a net force acting on the 
CME between 2-3 Rq that quickly rises to a peak total 
force of ~10 16 N and then falls rapidly. 

If we assume that the peak force observed for the 2008 
December 12 CME is the net force due to similar forces 
used in the above models, such as the solar wind drag, 
gravity, and some form of magnetic CME driver e.g., a 
J x B force, we may estimate their relative contribution. 
The force due to solar wind drag on the CME is given 

by 



F d 



CdPsw^-cme 



(v 



(2) 



where M cme is the CME mass, v is the CME velocity, 
Cd is the drag coefficient, p sw is the solar wind mass 
density, A cme is the CME area expos ed to solar wind 
drag and v s1 l , is the solar wind velocity ( |Maloney fc Gal- 
lagher^ 10 ). To estimate the effects of this force we use 
Psw = n p m p , where m p is proton mass, and assume ion- 



ization fraction of \ = 1 such that n p = n e [cm 



Elec- 



tron density, and hence proton density, is then given by 
an interplanetary density model derived fr om a special 
solut ion of the Parker solar wind equation ( |Mann et al.| 
1999), solar wind velocity values as a function ol height 
are also determined using this model. A cme is estimated 



using the expression derived in Byrne et al. ( |2010[ ) for 
latitudinal angular width of the CME as a function of 
height, A0i at (r) = 26r 22 . This is used to derive an arc 
length of the CME front and, as above, making the as- 
sumption A6i ong = 2x AOiat, the two arc lengths derived 
from these angles then give the surface that the solar 
wind acts on, thus A cme = 1352r 2 4 4 . Setting the drag 
coefficient Cd = 1, and using the Mann et al. (1999) 
model to derive a density and a solar wind velocity of 
2.3xl0 5 cm -3 and 70kms _1 , respectively, equation [1] 

then gives a force of Fd — —8.0 x 10 12 r N for solar wind 
drag at ~3i?0, where r is a unit vector in the positive 
radial direction. 

A simple estimate of force due to gravity is given by 
F g = GMQM cme /f 2 , where G is the universal gravita- 
tional constant, M is solar mass, is CME mass, 
and r is a heliocentric position vecto^J Given a CME 
mass of 3.4 xlO 15 g the force due to gravity at a heliocen- 
tric distance of 3 Rq is F g = -1.0 x 10 14 fN. The only 
remaining contribution is due to some form of magnetic 



4 Ideally the heliocentric distance of the CME centre of mass 
would be used here. However an unknown amount of mass is ob- 
scured by the coronagraphs occulting disk, making the mass dis- 
tribution and hence COM difficult to determine. Thus the CME 
front height is used in the calculation of force due to gravity 



8 



CME driver, F mag , which is estimated using 



-Fmag -^total 



i ~ F g (3) 

(the pressure gradient in the CME equation of motion 
is assumed to be negligable and has been omitted here). 
Using the above values, the total magnetic contribution 
to CME force is calculated to be F mag ~ 4.5xl0 14 fN at 
3 Rq , indicating that this is the largest driver of CMEs 
at low coronal heights. Lorentz force dominated dy nam- 



ics in early p hase CME propagation are reported in |Bein 
et al. ( 2Q11| ), in which a statistical study of a large sam- 
ple of CMEs in EUVI, COR1, and COR2 indicated an 
early phase acceleration for the majority of CMEs that 
is attributable to a Lorentz force. A similar result of 



an observational study by Vrsnak (2006) found that the 



Lorentz force plays a dominant role within a few solar 
radii. It should be noted that although we have labelled 
the force F magi there is no distinction on the exact form 
of this force e.g., whether it is magnetic pressure, mag- 
netic tension, or a Lorentz self-force that acts as the 
driver. Also, any no n-radial motion of t he CME, such 
as that described in Byrne et al. (2010), is not taken 



into account here; any force estimates are purely radial 
in direction. 

6. CONCLUSION 

The STEREO COR1/2 coronagraphs have been used 
to determine the mass development of the 2008 Decem- 
ber 12 CME. Knowledge of the longitudinal propaga- 
tion angle of the CME allowed for a significant reduc- 
tion in the mass uncertainty, giving a final estimate of 
3.4 zb 1 .0xl0 15 g. Using k inematics results of a previous 
study (|Byrne et al.||2010 ), the velocity and acceleration 
of the CME were combined with the mass measurements 



to determine the kinetic energy and total force on the 
CME. The early phase propagation of the CME was 
found to be dominated by a force of peak magnitude 
of 3.4±2.2xl0 14 N at ~3.Oi? , after which the magni- 
tude declines rapidly and settles to and average of 3.8 ± 
5.4xl0 13 N. This early rise and fall in total force (or ac- 
celeration) is in a greement with previous observations o f 



CME kinematics ( jGallagher et al.|2003||Bein et al.|20lT]) 
Simila rly re sults of observationa l studies by |Zhang et aT] 
( |2001[ ) and | Zhang et al.| (2004) also show early phase 
peak acceleration between 2-5 Rq and forces on the or- 
der of 10 15 N and 10 12 N. The kinetic energy shows an 



initial rise towards 6.3 ± 3.7 x 10 ergs at ~3 Rq, beyond 



which it rises steadily to 4.2 ± 2.5 xl(r u ergs at ~18i? , 
su ch order of mag nitudes are similar to those reported 
in |Vourlidas et al.| ( |2000| ); [ Emslie et al] (|2004[) and a re 
typical of CME kinetic energies ( [ Vourlidas eTaLpOlOl ) . 

Such CME kinematics and dynamics property esti- 
mates cannot be carried out when unknown propaga- 
tion angle hinders an accurate calculation of CME mass, 
hence adding unacceptable uncertainty to any subse- 
quent calculations. This highlights the need for similar 
studies using the STEREO mission's ability to accurately 
determine the physical properties of CMEs, such as mass, 
with remarkably reduced uncertainty. Increasing the ac- 
curacy of force estimates of other well studied CMEs will 
allow for a more complete view of the magnitude of the 
forces influencing CME propagation and will allow model 
parameters to be more accurately constrained. 

This work is supported by the Irish Research Coun- 
cil for Science, Engineering and Technology (IRC- 
SET). We also extend thanks and appreciation to the 
STEREO /SECCHI consortium for providing open access 
to their data. 
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